Suppression of Ground-State Magnetization in Finite-Sized Systems Due to 

Off-Diagonal Interaction Fluctuations 
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We study a generic model of interacting fermions in a finite-sized disordered system. We show that 
the off-diagonal interaction matrix elements induce density of states fluctuations which generically 
favor a minimum spin ground state at large interaction amplitude, U. This effect competes with 
the exchange effect which favors large magnetization at large U, and it suppresses this exchange 
magnetization in a large parameter range. When off-diagonal fluctuations dominate, the model 
predicts a spin gap which is larger for odd-spin ground states as for even-spin, suggesting a simple 
experimental signature of this off-diagonal effect in Coulomb blockade transport measurements. 
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Ferromagnetic instabilities result from the combined 
effect of the electronic interactions together with the 
Pauli principle. The interaction energy can be minimized 
when the fermionic antisymmetry requirement is satisfied 
by the spatial wavefunction leading to alignment of spins 
and a large ground-state spin magnetization (a familiar 
example of this is Hund's rule in atoms). In contrast 
when the interaction is weak minimal spin is favored be- 
cause it costs kinetic energy to flip a spin as it must then 
be promoted to a higher energy level. When treating 
ferromagnetism in metals, because of the locality of the 
Pauli principle, magnetic instabilities are usually studied 
in the framework of the Hubbard model taking only the 
short-range part of electronic interactions into account, 
so that only pairs of electrons of opposite spin interact. 
As the magnetization increases the number of interacting 
pairs decreases and the system spontaneously magnetizes 
for sufficiently strong interactions. In a finite-sized sys- 
tem this Stoner instability (l) occurs when the typical ex- 
change interaction between two states close to the Fermi 
energy is equal to the one-particle level spacing which for 
a clean system with Hubbard interaction gives U c = A. 

There has been much recent interest in the Stoner in- 
stability of finite-size disordered metals such as quantum 
dots and metallic nanoparticles Building on earlier 
perturbative work || Andreev and Kamenev Q recently 
found a significant reduction of the Stoner threshold in 
disordered systems due to spatial correlations in diffusive 
wavefunctions which enhance the average exchange term. 
More recently, Brouwer et. al. (5| considered the effect of 
mesoscopic wavefunction fluctuations and found an asso- 
ciated increase of the probability of non-zero ground state 
spin magnetization below the Stoner threshold. Within 
a similar model, Baranger et. al. || proposed that spon- 
tanous magnetization effects could explain kinks in the 
field-dependence of Coulomb blockade resonances. The 
purpose of the present paper is to point out a compet- 
ing effect of interactions which suppresses the probability 
of ground-state magnetization and has not been treated 
in any of the previous works on itinerant magnetism 
of disordered systems. The mean-field treatments lead- 



ing to an exchange term in the effective hamiltonian for 
disordered metals neglects the effects of off- diagonal in- 
teraction matrix elements [Q; however it is well-known 
from studies of nuclei and atoms (?]] that the band-width 
of the many-body density of states in finite interacting 
fermi systems is actually determined by the fluctuations 
of these off-diagonal matrix elements. When one intro- 
duces the spin degree of freedom into these models we 
shall see below that one immediately finds that these 
fluctuations are largest for the states of minimal spin. 
This effect then strongly increases the probability that 
the extremal (low-lying) states in the band are those of 
minimal spin and opposes the exchange effect. We ex- 
pect this effect to be significant in quantum dots and to 
suppress the possibility of high spin ground states. 
We start from the Hamiltonian for n spin- 1/2 particles 



(1) 



s' ' =t,l are spin indices. The m/2 different one-body 
energies are distributed as e a 6 [—m/2; m/2] so as to fix 
A = 1 with spin degeneracy. The interaction commutes 
with the z-component a z of the total magnetization a so 
that the Hamiltonian acquires a block structure where 
blocks are labelled by <j z and due to Spin Rotational Sym- 
metry (SRS) subblocks of given a > \<r z \ appear within 
each of these blocks. Each block's size is given in term 

of binomial coefficients as N(a z ) — J )■ 

The Hamiltonian (1) can be viewed as a generic model 
of interacting fermions expressed in the basis of Slater 
determinants constructed from the eigenstates ip a of the 

t 



corresponding free fermions model Hq = ^2e a c' a s ^ a ^ s . 
In this basis the interaction matrix elements arc given by 

J a,B 



UZ'" n = fdrdr'U(r-r')ipZ(^p(r')<ip y {r)-ips(r'), where 



U(f — r') is the interaction potential. Due to disor- 
der or chaotic boundary scattering the wavefunctions ip a 
have a random character leading to fluctuations in U^'g 
around their average value. We take these fluctuations 
to be random with a zero-centered gaussian distribution 
of width U. This gives a contribution H similar to the 
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second term in the right-hand side of (|l|) with a distribu- 
tion P(U^'t) ex e~^ Ua -^ ^ 2U of interaction matrix ele- 
ments. Only diagonal matrix elements U^'p and U^'p 
have a nonzero average leading to mean-field charge- 
charge and spin-spin diagonal interactions Q . We neglect 
the charge-charge contribution as it has no influence on 
the magnetization and this leaves us with the following 
effective Hamiltonian 

n = H a + H-\UJ2^sp (2) 

s a = J2 S t c a s ^s,tC a ,t are spin operators and the ferro- 
magnetic spin-spin interaction has a strength XU > 0. 
Without it (A = 0), the Hamiltonian (0) within each 
spin block is precisely the Two-Body Random Interac- 
tion Model (TBRIM) introduced in nuclear physics ||, 
and used to study thermalization || and the emergence 
of quantum chaos in few-body systems and statisti- 
cal features observed in shell model calculations ||lT| . |l2| . 
A similar model has been shown recently to be consistent 
with the observed gaussian distribution of peak-spacings 
in Coulomb blockade resonances through quantum dots 
|U|. One key feature of the TBRIM is that the many- 
body density of states (MBDOS) has an approximately 
gaussian shape with a variance proportional to the con- 
nectivity K, i.e. the number of non-zero matrix elements 
in each row 0. Similarly we can estimate the variance 
of the MBDOS of the Hamiltonian (§) for fixed (a,a z ) 
and U/A > 1 as 

-J— J2 HIj8{*W _ Uz)5{u {D _ a) „ KU 2 (3) 
N(a z ) 

where Sj j — (I\H\J) and \I) refers to a Slater determi- 
nant. Hence each block's bandwidth goes as \f~KU with 
a CT z -dependent connectivity which can be expressed as 
K(n,m,a z ) w l + C(n/2 + a z ,m/2) + C(n/2~a z ,m/2) + 
l/2((n/2) 2 -a 2 )((m/2-n/2) 2 -a 2 ) in term of the func- 
tion C(n, m) — n{m — n) + n(ri — l)(m— n)(m — n— l)/4. 
The factor 1/2 in front of the last contribution to K is 
needed to take into account the effect of SRS. In the di- 
lute limit approximately half of the spin-flip transitions 
which conserve a z are not allowed because they don't 
conserve total spin (e.g. change a singlet to triplet). The 
estimate (Jfy assumes that each matrix element has the 
same variance, which for a generic off-diagonal element is 
~ U 2 . However it is easily seen that diagonal matrix ele- 
ments Hi j have an enhanced variance ~ (3n 2 / 1 A+a 2 ) 2 U 2 
which induces deviations from (H) for large filling and 
large magnetization. Nevertheless, these matrix elements 
can be neglected in the dilute and weakly-magnetized 
limit 1 <C n/2 + <j z <C m where the larger number of off- 
diagonal matrix elements dominates the variance when 
(m — n) 2 3> (3n 2 /4 + a 2 ). In the right inset to Fig. 1 
we show plots of K(n,m,o~ z ) for different filling factors 
v = n/m as well as a comparison with the true variance 



of the MBDOS for v = 3/8. Deviations from the esti- 
mate (^) are small, even at this rather large filling and 
increase with increasing magnetization in agreement with 
the above reasoning. Moreover in the left inset to Fig. 1 
we show that the full MBDOS follows the scaling (||) with 
significant deviations only at large magnetization, so that 
it is plausible that this scaling also determines the tails 
in which the ground state energies will be found. 




E/U 

FIG. 1. Density of states for the Hamiltonian H with n — 6 
particles and m = 16 orbitals, corresponding to the magneti- 
zation blocks a z — —3 (solid line), -2 (dotted line), -1 (dashed 
line) and (dotted-dashed line). Left inset: rescaled density 
of states showing the approximate scaling in E/K 1/I2 U. Right 
inset : Normalized connectivity K(a z )/ K(0) vs. magnetiza- 
tion for filling factors v = 1/10 (solid line), 3/8 (dashed line) 
and 1/2 (dotted-dashed line). We compare this estimate to 
the true variance for v = 3/8 obtained numerically from the 
data of the main figure (diamonds). 

Our essential finding follows from the simple features 
of the model already stated. The full MBDOS is a sum 
of approximately gaussian contributions from each spin 
block with a variance proportional to the corresponding 
connectivity. The latter is a monotonously decreasing 
function of a z . Hence the broadest MBDOS corresponds 
to the minimally magnetized block and the ground state 
will be found in this block with increased probability . 
Assuming, as just discussed, that the tails of the distribu- 
tion scale with the variance with a factor (3 and neglecting 
contributions arising from Hq 7 the typical spin gap can 
be estimated (for A = 0) as 

I) - yjK {\u min \ + 1)] (4) 

This multiple gaussian structure of the MBDOS and the 
scaling with \f~K obtained from numerical calculations 
are shown in Fig. 1 for A = 0. For A ^ 0, the spin-spin 
interaction induces relative shifts of each block's MBDOS 
which eventually will shift the finite spin blocks suffi- 
ciently to overcome the larger fluctuations of the minimal 
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spin MBDOS. This is the competition between exchange 
and off-diagonal fluctuations already mentioned. How- 
ever for reasonable values of A the off-diagonal fluctua- 
tions strongly reduce the probability of exchange- induced 
magnetization (see Fig. 3). 
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FIG. 2. Dependence of the finite-size spin gap in the num- 
ber n of particles. Points correspond to numerical results for 
m = 10 (full circles), 12 (empty squares), 14 (full diamonds) 
and 16 (empty triangles) and the solid line to the dilute es- 
timate (0) with a numerical factor /3 = 1.5. For the case 
m — 16 and 1000 Hamiltonian realizations, the error bars 
indicate the r.m.s. of the gap distribution while the dashed 
and dotted-dashed lines show the numerically computed vari- 
ances (Left-hand side of eq. (3)) for the full Hamiltonian and 
after setting to zero non-generic interaction matrix elements 
respectively. 

In Fig. 2 we show the computed spin gap between 
the minimally magnetized ground-state and the first spin 
excited level for A = in the limit of dominant interac- 
tion, i.e. neglecting Hq in (|l|). One of the main features is 
a strong even-odd effect which is reminiscent of a similar 
behavior in the limit of vanishing interactions. However 
the origin here is the fluctuating interaction and the en- 
ergy differences scale as U instead of A. We next note 
that the gap first increases with increasing number of par- 
ticles before it seems to stabilize above n = 6. We have 
checked (dashed and dot-dashed lined in Fig. 2) that this 
behaviour, which is not captured by the dilute estimate 
(Q), is partly due to the neglect in (j|) of nongeneric ma- 
trix elements with enhanced variance mentioned above. 
However, even though the exact variance gives a much 
better estimate, it still underestimates the gap at larger 
n and we have numerically determined that this is due 
to a strong positive correlation of the ground state ener- 
gies in adjoining spin blocks (~ 0.9). Such correlations, 
although interesting, are not suprising since the differ- 
ent block hamiltonians are not statistically independent 
(many of the same two-body matrix elements appear in 
both). 



We next switch on the mean-field spin-spin interaction 
A > which induces energy shifts of —XU \<r z \ (\a z \ + 1) 
0. On average the spin gap becomes A s = A^ — XU, 
where A = (5— (— l)")A/2, i.e. the relative shift between 
the two lowest magnetized blocks is larger for odd number 
of particles. The variance of the gap distribution is unaf- 
fected and we can already conclude that the probability 
P(a z > 0) of finding a magnetized ground-state is re- 
duced by the off-diagonal matrix elements and saturates 
above a finite value U c , since the width of the gap distri- 
bution is proportional to its average ~ U. This is shown 
on Fig. 3 where P(cr z > 0) is plotted against U/A for 
different values of A. On the same graph we show numer- 
ical results obtained after setting to zero the off-diagonal 
matrix elements. The data unambiguously reflect the 
strong demagnetizing influence of the off-diagonal matrix 
elements. 
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FIG. 3. Probability for a magnetized ground-state as a 
function of U/A for 2000 realizations of Hamiltonian (Q) with 
n = 5 and m — 10, A = 2, 4, 5 and 7 (full symbols, from bot- 
tom to top). Empty symbols show the corresponding curves 
after setting to zero the off-diagonal matrix elements. Stars 
correspond to n — 5, m — 14 and A = 7. 

The physically relevant value of A will depend on 
the microscopic details of the system. Indeed A is 
given by half the ratio of the mean exchange interaction 
with the fluctuations of off-diagonal matrix elements Q] 
A= (U^p)/2 r.m.s. Semiconductor quantum dots 
with poor screening and extended, chaotic single-particle 
wavefunctions should have Awl while for extended dif- 
fusive metallic systems, we get a much larger spin-spin 
interaction A ~ g where g > 1 is the system's conduc- 
tance ||. From Fig. 3 however, we see that even in 
the regime of dominating interactions, a nonmagnetized 
ground state is more probable for A < 5 {XU <C A^). 

We finally consider the influence of an external mag- 
netic field which only introduces a Zeemann coupling. 
This situation can be experimentally realized by applying 
a magnetic field in the plane of a two-dimensional elec- 
tron gas. This Zeemann term does not affect the Hamil- 
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tonian's block structure, but only shifts each block's MB- 
DOS by an amount gnsBa z . Due to the spin gap dis- 
cussed above, a finite magnetic field of average magnitude 
(B c ) = A s /g[iB is necessary to magnetize the system. 
The even-odd effect emphasized in Fig. 2 results in a 
critical field to flip one spin which is significantly larger 
for a a z — ±1/2 (odd) ground state as compared to a 
a z = (even) ground state. More generally the gap for 
a lower spin state is smaller than that for the state of 
higher spin (this follows from the right inset to Fig. 1). 
This aspect of our theory can be tested experimentally 
by studying the in-plane magnetic field dependence of 
the position of Coulomb blockade conductance peaks at 
very low temperature r« A. The resonant gate voltage 
is given by a difference of two many-body ground-state 
energies eV g n = E® +1 — E®, and it is always the difference 
of an even-odd pair. The peak position behaves like 



eV Q n (B) = E° n 



E°+ gf i B BSa z (n) 



(5) 



where Sa z (n) is the magnetization difference between the 
two consecutive ground-states. Without magnetization 
6<j z (n) = (-l) n /2 and one has \dV g /dB\ = g/j, B /2. As 
B is increased the ground-state with even number of elec- 
trons is most likely to magnetize first, exactly reversing 
the slope of two consecutive peaks; then as the field in- 
creases further the odd state will likely flip, restoring the 
original slope. As long as consecutive ground states never 
differ by more than one unit of spin the absolute value of 
the slope will remain constant as the system polarizes. 
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FIG. 4. Schematic of the conductance peaks in a 2D quan- 
tum dot as a function of an in-plane magnetic field for m — 14, 
A = 1 (left) and 2.5 (right) and U/A = 4 corresponding to 
the addition of the n — 3,4, ...10* electrons (from bottom to 
top). Dotted-dashed lines indicate slopes of ±g/i_g/2. Larger 
slopes for which spin-blockade effects strongly reduce the peak 
height [15] are indicated by dashed segments. Note that due 
to the subtraction of the average charge-charge interaction, 
the model does not reproduce the charging energy so that the 
vertical distance between consecutive peaks is arbitrary. 

However if there exist many magnetized ground states 
then one expects a range of slopes to occur. In this case 
the corresponding peak heights will be strongly reduced 



by the spin blockade mechanism (see dashed lines on Fig. 
4) |l5| which should be easily visible experimentally. This 
argument neglects changes in the g-factor of the electron 
with changing n, which presumably are slow. The even- 
odd behavior of A s is qualitatively similar to the non- 
interacting case, however the scale in gfJ-sB over which 
spin flips occur is determined by U and not A. Typically, 
this results in an increase of the field necessary to achieve 
full polarization on the dot. This is illustrated on Fig. 4 
where the peaks positions are drawn as a function of the 
Zeemann coupling for A = 1.5 and 3. It is clearly seen 
that at small A, \dV g /dB\ is constant and corresponds to 
a minimal 6<j z , while increasing A gives different slopes 
in agreement with the above reasoning. 
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